par(mfrow=c(3, 2))
data1_8<-read.csv("data/时间序列分析——基于R（第2版）案例数据/csv/A1_8.csv")
ts_overshort<- ts(data1_8$overshort)

# 时序图
plot(ts_overshort, type="o", pch=5, col='#39CBB4')

# 自相关图
acf(ts_overshort, lag.max=30)
pacf(ts_overshort, lag.max=30)

# ADF检验
install.packages("aTSA")
library(aTSA)
adf.test(ts_overshort, nlag=3)

# 白噪声，纯随机检验
for( k in 1:3) print(Box.test(ts_overshort, lag=6*k, type="Ljung-Box"))

# 